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Executive  Summary 

Metal-oxide-semiconductor  (MOS)  structures  There  is  an  active  controversy  over  whether 
are  being  used  with  increasing  frequency  in  inte-  avalanching  is  important  in  oxide  breakdown.  The 


grated  circuits  for  military  systems.  They  will  be  lo¬ 
cated  in  systems  where  they  may  be  exposed  to 
electrical  overstress  from  the  electromagnetic 
pulse  (EMP)  and  high  power  microwaves  (HPM’s). 
Because  of  the  thinness  of  the  oxide  layer,  it  is  pos¬ 
sible  that  damage  may  occur  in  the  oxide  at  lower 
levels  than  in  active  semiconductor  junctions. 

Oxide  breakdown  is  becoming  the  dominant 
damage  mechanism  in  integrated  circuits  as  the 
oxide  thicknesses  are  decreased  in  step  with  de¬ 
vice  dimensions.  The  Army’s  concern  with  this  type 
of  damage  is  most  relevant  in  the  very-high-speed 
integrated  circuit  (VHSIC)  program.  No  known  com¬ 
parable  studies  are  being  conducted  on  this 
problem. 


Harry  Diamond  Laboratories  DIODE  computer  pro¬ 
gram  has  been  used  to  study  avalanche  breakdown 
in  silicon  dioxide.  The  calculations  strongly  support 
avalanching  as  the  dominant  cause  of  breakdown, 
at  least  for  thicknesses  greater  than  50  nm.  The 
calculations  are  limited  to  maximum  breakdown 
times  of  the  order  of  100  ps,  whereas  the  only  pub¬ 
lished  experimental  measurements  are  in  the  1-  to 
100-s  time  scale,  or  breakdown  by  laser  irradiation. 

There  are  no  present  plans  to  continue  the 
study  of  avalanche  breakdown  in  Si02.  However,  it 
is  strongly  urged  that  an  experimental  study  be 
made  for  Si02  breakdown  at  rates  of  voltage  rise 
comparable  to  those  used  in  the  calculations. 


1.  Introduction 

Breakdown  of  silicon  dioxide  films  has  become 
an  increasing  problem  as  the  films  have  decreased 
in  thickness  in  parallel  with  the  device  dimensions 
in  very-large-scale  integration.  Although  break¬ 
down  in  Si02  has  been  studied  for  several  decades, 
there  is  still  an  active  controversy  over  the  physical 
mechanism  of  that  breakdown. 

Many  early  investigators  proposed  an  ava¬ 
lanche  (impact  ionization)  breakdown  mechanism, 
as  proven  in  gases  and  later  in  semiconductors. 
O’Dwyer  [1]*  proposed  that  the  avalanche  pro¬ 
cess  was  augmented  by  the  space  charge  effects 
of  slow-moving  holes.  DiStefano  and  Shatzkes  [2] 
obtained  agreement  with  the  measured  depend¬ 
ence  of  the  breakdown  field  upon  oxide  thickness 
with  an  avalanche-recombination  model.  However, 
this  model  uses  the  ionization  band  gap  and 
electron-phonon  scattering  length  and  is  not 
suitable  for  routine  calculations.  An  extensive  list 
of  references  is  presented  in  this  report  [2], 

Another  physical  process  suggested  as  impor¬ 
tant  in  oxide  breakdown  is  tunneling — either  direct 
tunneling  or  Fowler-Nordheim  tunneling  [3],  Chen 
et  al  [4]  agree  with  measured  breakdown  times 
(ranging  approximately  from  1  x  lo2to3  x  104  s)as 
a  function  of  applied  field  by  assuming  75-percent 
tunneling  and  a  25-percent  avalanche  contribution. 

Weinburg  [5]  used  a  carrier  separation  method 
which  he  interprets  as  evidence  against  impact 
ionization.  Recently,  DiMaria  et  al  [6]  used  this 
technique  and  two  others  to  reach  the  same  con¬ 
clusion  for  oxide  films  of  less  than  7  nm  (70  A). 

Although  still  controversial,  the  evidence  for 
avalanche  breakdown  is  sufficient  to  justify  an  in¬ 
vestigation  of  oxide  breakdown  with  the  use  of  the 
Harry  Diamond  Laboratories  (HDL)  DIODE  com¬ 
puter  program.  Section  2  discusses  the  input 
parameters  used  in  the  calculations,  and  the 
following  sections  give  results  of  calculations  for  a 
constant  applied  voltage,  for  constant  voltage  ramp 
rates,  and  for  an  applied  sinusoidal  voltage.  Sec¬ 
tion  6  includes  further  discussion  and  conclusions. 


2.  Computer  Program  and  Input  Parameters 

The  HDL  DIODE  program  was  originally  used  to 
study  electrical  breakdown  in  gases  and  has  been 
used  extensively  to  study  second  breakdown  in 
silicon  [7],  The  voltage  is  applied  to  the  oxide  film 
through  an  RC  circuit.  For  these  calculations  the 
series  resistance  was  2  ohms  and  the  shunt 
capacitance  was  1  pF,  giving  a  voltage  risetime  of 
2  ps.  The  applied  voltage  was  an  initial  step  or  a 
portion  of  a  sinusoidal  voltage.  The  latter  option 
allows  a  constant  voltage  ramp  to  be  applied. 

The  most  important  calculation  parameters  for 
avalanche  breakdown  are  the  ionization  coeffi¬ 
cients  for  the  electrons  and  holes.  In  Si02,  the  ef¬ 
fective  hole  mobility  is  a  factor  of  about  1  x  io6 
lower  than  that  for  electrons,  due  to  hole  trapping 
and  probably  a  hopping  process;  therefore,  ioniza¬ 
tion  by  holes  is  improbable.  Unfortunately,  there  is 
not  even  order-of-magnitude  agreement  on  the 
dependence  of  the  electron  ionization  coefficient 
upon  electric  field.  The  most  comprehensive 
survey  of  the  high  field  electronic  properties  of 
Si02  has  been  made  by  Hughes  [8].  Hughes  com¬ 
pares  the  ionization  coefficients  deduced  from 
laser  breakdown  experiments  in  bulk  material  and 
those  determined  from  near-dc  measurements  in 
thin  films.  There  is  rough  agreement  at  fields  about 
1  x  io7  Wcm,  but  orders-of-magnitude  disagree¬ 
ment  for  fields  of  half  that  value.  In  the  calculation 
of  the  ionization  coefficients  from  the  laser 
breakdown,  it  was  assumed  that  the  electron  densi¬ 
ty  increased  from  108cm-3  (one  electron  in  the 
volume  irradiated)  to  the  plasma  density  of 
1018cm-3.  Obviously  both  of  these  numbers  are 
only  approximate,  but  because  the  ionization  coef¬ 
ficient  appears  in  an  exponential,  its  error  is 
markedly  less.  Some  laser  breakdown  data  from  a 
table  in  Hughes'  paper  [8]  are  reproduced  in  table 
1.  It  is  seen  that  a  higher  field  is  required  for 

Table  1.  Measured  Breakdown  In  Si02 


Pulse 

Film 

Peak 

Method 

width 

(ns) 

thickness 

(nm) 

field 

(MV/cm) 

Reference 

Laser 

0.03 

— 

17 

9 

Laser 

5 

— 

7.4 

10 

breakdown  with  a  30-ps  pulse  width  [9]  than  is  re-  106  cm-3  at  300  K.  Calculations  were  made  for  an 
quired  for  a  5-ns  pulse  width  [10].  These  meas-  area  of  10_5cm2. 
urements  were  made  by  different  investigators,  but 

the  same  wavelength  of  1  nm  was  used.  3.  Constant  Applied  Voltages 


The  scatter  of  measured  breakdown  fields  in 
thin  films  is  roughly  ±20  percent  at  100  nm  and 
greater  for  thinner  films.  Also  shown  in  table  1  is  a 
thin  film  measurement  by  Harari  [11]  and  a  calcula¬ 
ted  result  [2]  that  best  fit  the  experimental  data. 
These  two  breakdown  fields  were  chosen  to  derive 
ionization  coefficients  to  use  in  the  calculations  of 
this  report.  The  form  of  the  ionization  coefficient,  a, 
was  assumed  to  be 


a  -  A  exp  (-B/E) , 


(D 


where  E  is  the  electric  field  and  A  and  B  are  the 
parameters  to  be  fitted.  A  further  assumption  was 
that  ad  =  1  for  breakdown,  where  d  is  the  oxide 
thickness.  This  assumption  is  true  only  for  equal 
electron  and  hole  coefficients  and  is  often  made 
even  when  undoubtedly  not  true.  The  parameters 
were  calculated  to  be  A  *  1  x  io7  cm-1  and  B  « 
4.6  x  107  cm/V.  Chen  [4]  used  A  -  3.3  x  106  cm-1 
and  B  *  7.8  x  107  cm/V.  Assuming  an  electron 
saturation  velocity  of  2x10?  cm/s  as  given  by 
Hughes  [8],  the  ionizations-per-second  data  from 
the  laser  measurements  can  be  converted  to  the 
ionizations-per-centimeter  form  of  equation  (1). 
This  gives  A  -  1  x  io6cm-1  and  B»  4.6  x  io7  cm/V, 
precisely  one  order  of  magnitude  smaller  than 
those  used  in  this  report.  The  data  used  by  Chen  [4] 
agree  with  the  laser  data  at  a  field  of  approxi¬ 
mately  2  x  io7  V/cm,  but  are  an  order  of  magnitude 
lower  at  a  field  just  under  1  x  io7  V/cm.  An  order  of 
magnitude  change  in  the  ionization  coefficient 
typically  makes  only  a  25-percent  change  in 
breakdown  field,  a  typical  scatter  in  published 
measurements. 

Second  in  importance  in  avalanche  breakdown 
calculations  are  the  electron  and  hole  mobilities 
and  saturation  velocities.  Hughes  [8]  gives  a  low 
field  mobility  of  electrons  in  SiO,  of  21  cm2/V-s,  and 
a  saturation  velocity  of  2x  10'cm/s.  He  gives  a 
constant  hole  mobility  of  2  x  icr5  cm2/V-s  up  to  the 
highest  field  measurement  (6xio6V/cm).  These 
values  were  used  in  the  present  calculations.  Other 
parameters  used  in  the  calculations  were  a  dielec¬ 
tric  constant  of  3.9  and  an  intrinsic  density  of 


Computer  calculations  were  made  simulating 
breakdown  in  a  100-nm  (1000-A)  oxide  layer  with 
constant  applied  voltages.  Figure  1  shows  the 
results,  where  the  current  density  is  plotted  as  a 
function  of  time.  The  curves  of  this  figure  are 
qualitatively  similar  to  those  measured  in  semicon¬ 
ductors  and  gases.  The  dc  breakdown  voltage  is 
100  V  for  this  oxide  thickness,  as  determined  from 
the  ionization  coefficients.  Thus,  overvoltages  of 
145  to  180  percent  were  needed  to  obtain  break¬ 
down  in  the  100-ps  time  frame.  Since  the  electron 
transit  time  across  this  oxide  film  is  0.5  ps,  it 
becomes  expensive  to  extend  calculations  to  much 
greater  times. 

There  are  no  published  experimental  data  which 
can  be  compared  directly  with  figure  1.  However, 
Harari  [11]  has  measured  current-voltage  charac¬ 
teristics  in  thin  oxide  films  with  slow  rising  voltage 
ramps.  Selected  curves  (three  of  five  oxide 
thicknesses)  from  one  of  his  figures  are  repro¬ 
duced  in  figure  2.  In  figure  2(b),  the  voltage  was 
ramped  to  a  given  final  voltage,  at  which  it  was 
held.  The  current  decayed  at  the  constant  voltage 
until  breakdown  occurred.  Harari  postulated  that 


Figure  1.  Computer  current  density  growth  as  s  function  of 
time.  Parameter  Is  applied  voltage. 
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VOLTAGE (V) 

Figure  2.  Selected  experimental  current  voltage  curves 
from  Hararl  [11].  (a)  Constant  ramp  rates:  (1)  0.013  V/s, 
(2)  0.15  V/s,  (3)  1.69  V/s,  (4)  30  V/s.  (b)  Ramp,  followed  by 
constant  voltage. 

the  current  decay  was  due  to  electron  trapping.  An 
alternative  explanation  is  that  the  decay  was  due  to 
electron-hole  recombination. 

The  curves  in  figure  1  were  calculated  with  no 
recombination,  i.e.,  an  infinite  lifetime.  Further 
calculations  were  made  with  various  recombina¬ 
tion  lifetimes.  The  results  are  shown  in  figure  3.  The 
applied  voltage  is  260  V.  The  slight  decrease  in 
breakdown  voltage,  as  compared  to  figure  1, 
resulted  from  using  only  25  spatial  grid  points  in¬ 
stead  of  50.  This  reduced  the  computing  cost  by  a 
factor  of  four.  It  is  seen  in  figure  3  that  increasing 
the  recombination  rate  (reducing  the  lifetime) 
decreases  the  current  growth  rate.  The  sharp  cur¬ 
rent  drop  seen  in  the  250-ms  lifetime  curve  is  due 
to  external-circuit-produced  relaxation  oscillations 
[7].  The  100-ms  lifetime  curve  shows  a  slight  cur¬ 
rent  decay  not  associated  with  the  external  circuit. 
This  shows  that  recombination  of  electrons  and 
holes  is  a  plausible  explanation  for  the  current 
decays  noted  by  Harari. 


TIME  (ps) 

Figure  3.  Calculated  current  density  growth  with  recombi¬ 
nation  lifetime  given  by  labeled  parameter.  Applied  voltage 
Is  260  V. 


The  lifetimes  used  in  figure  3  are  those  ap¬ 
plicable  at  low  current  densities,  i.e.,  for  carrier 
densities  near  the  intrinsic  density  of  106cm~3. 
The  effective  lifetime  varies  inversely  with  the  car¬ 
rier  density.  For  these  calculations  hole  densities 
exceed  1019cm~3  at  breakdown  (maximum  elec¬ 
tron  densities  are  strongly  dependent  upon  the  life¬ 
time  used).  The  extremely  short  effective  lifetimes 
explain  the  current  decay  in  the  picosecond  time 
scale  when  millisecond  intrinsic  lifetimes  are  used. 

The  low  mobility  of  holes  means  that  they  are 
essentially  stationary  in  the  time  scale  of  these  cal¬ 
culations,  whereas  many  electrons  exit  the  oxide  at 
the  anode.  This  causes  a  buildup  of  the  hole  space 
charge.  Calculated  electric  fields  for  the  125-ms 
lifetime  calculation  of  figure  3  are  shown  in  figure 
4.  The  negative  slope  of  the  electric  field  is  propor¬ 
tional  to  the  excess  hole  density.  The  holes  are 
created  closer  to  the  cathode  as  the  cathode  field 
increases  and  the  anode  field  drops.  A  small  in¬ 
crease  in  the  field  means  a  large  increase  in  the 
ionization  rate,  since  the  ionization  coefficient  in¬ 
creases  exponentially  with  the  field. 
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Figure  4.  Field  distributions  as  a  function  of  time  for  125-ms 
lifetime  curve  of  figure  3. 


Avalanche  breakdown  at  moderate  overvolt¬ 
ages  requires  some  feedback  mechanism.  In  sili¬ 
con,  ionization  by  holes  as  well  as  space  charge  ef¬ 
fects  causes  the  feedback.  In  Si02  there  is  essen¬ 
tially  no  hole  ionization  so  the  field  distortion  pro¬ 
vides  the  only  feedback  mechanism.  For  low  re¬ 
combination  rates  and  high  overvoltages,  the  hole 
space  charge  is  confined  to  the  anode  region.  For 
higher  recombination  rates  and  lower  overvolt¬ 
ages.  the  hole  space  charge  maximum  moves  to 
the  cathode.  This  may  explain  why  some  exper¬ 
imenters  report  finding  an  accumulation  of  holes  at 
the  cathode,  whereas  others  do  not. 

Calculations  with  constant  applied  voltages 
were  also  made  for  50-  and  200-nm  oxide  thick¬ 
nesses.  The  results  are  shown  in  figure  5,  where 
the  current  is  plotted  as  a  function  of  the  oxide  volt¬ 
age.  The  solid  curves  were  calculated  with  a  re¬ 
combination  lifetime  in  intrinsic  oxide  of  500  ms. 
The  260-V  dashed  curve  shows  a  current  decay 
before  breakdown,  the  latter  defined  by  the  appear¬ 
ance  of  a  negative  differential  resistance.  These 
curves  are  qualitatively  similar  to  those  measured 
by  Harari  (fig.  2)  even  though  there  is  a  difference 
in  time  scale  of  12  orders  of  magnitude.  The  nomin¬ 
al  dc  breakdown  voltages,  determined  by  the  ion¬ 
ization  coefficients  and  the  criterion  that  ad  -  1, 


Figure  5.  Current  density  es  a  function  of  voltage  for  three 
oxide  thicknesses.  Parameter  for  each  curve  Is  applied 
voltage.  Recombination  lifetime  is  500  ms  for  solid  curves 
and  100  ms  for  dashed  curves. 


are  59,  100,  and  174  V  for  the  oxide  thicknesses  of 
50,  100,  and  200  nm,  respectively.  The  cost  of  cal¬ 
culations  is  inversely  proportional  to  the  oxide 
thickness,  so  thinner  oxide  calculations  have  not 
been  made.  The  ad  =  1  criterion  for  breakdown 
precludes  breakdown  for  d  less  than  10~7  cm 
(1  nm)  with  the  ionization  coefficients  used  herein. 

4.  Constant  Voltage  Ramps 

Calculations  have  also  been  made  with  an  ap¬ 
plied  voltage  increasing  linearly  with  time,  a  cons¬ 
tant  voltage  ramp  rate.  The  voltage  ramp  is  pro¬ 
duced  by  the  initial  portion  of  a  sinusoidal  applied 
voltage.  An  initial  approach  voltage  was  used  to  re¬ 
duce  computer  cost.  The  initial  voltage  was  varied 
to  ascertain  that  the  calculated  breakdown  (max¬ 
imum)  voltage  was  unchanged.  The  calculations 
shown  in  figure  6  were  made  with  a  dc  voltage  of 
200  V  added  to  the  sinusoidal  voltage.  The  initial 
oxide  voltage  was  180  V.  It  is  seen  that  the  break¬ 
down  voltage  decreases  as  the  ramp  rate  de¬ 
creases.  Calculations  in  silicon  [7]  have  shown  that 
the  overvoltage  (ramp  voltage  maximum  minus  the 
dc  breakdown  voltage)  increases  with  the  2/3 
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Figure  6.  Voltage  as  a  function  of  time  with  various  ramp 
rates  as  given.  Initial  voltage  was  180  V. 


power  of  the  ramp  rate.  Accordingly,  trial  dc  break¬ 
down  voltages  were  chosen  to  plot  overvoltage  ver¬ 
sus  ramp  rate  for  the  data  of  figure  6.  The  results 
are  shown  in  figure  7(a).  A  dc  breakdown  voltage  of 
230  V  gives  the  best  straight  line  fit.  The  slope  of 
the  straight  line  is  0.22.  Calculations  were  also 
made  for  50-  and  200-nm  oxide  thicknesses,  where 
the  straight  line  slopes  were  0.19  and  0.28,  respec¬ 
tively.  For  the  200-nm  oxide  thickness,  calculations 
were  also  made  for  various  recombination 
lifetimes.  The  dc  breakdown  voltages  needed  to  ob¬ 
tain  straight  lines  increased  as  expected  as  the  life¬ 
time  decreased,  but  no  definite  change  of  slope 
was  found. 

The  experimental  data  of  Harari  (fig.  2(a))  were 
used  to  plot  the  overvoltage  against  the  ramp  rate 
in  figure  7(b).  Straight  lines  with  slopes  of  0.23  to 
0.33  were  found  with  appropriately  assumed  dc 
breakdown  voltages.  There  also  seems  to  be  a  sys¬ 
tematic  variation  of  the  intercepts  of  the  curves 
with  oxide  thickness  for  both  the  calculations  and 
the  measurements.  The  extrapolation  of  the  50-nm 
straight  line  to  zero  on  the  log  dV/dt  scale  (almost 
12  orders  of  magnitude)  gives  an  overvoltage  of 
about  0.4  V.  This  is  not  incompatible  with  Harari’s 
data  as  shown  in  figure  8.  The  dc  bref  'xlown  volt¬ 
ages  obtained  from  the  straight  line  fit  of  the  over- 
voltage-ramp  rate  plots  of  Harari’s  data  reduce  the 
breakdown  voltage  for  a  10-nm  oxide  from  20  to 
16  V.  These  reduced  voltages  give  closer  agree¬ 
ment  with  the  breakdown  voltages  found  for  thicker 
oxides  by  DiStefano  [2]. 


(a) 


RAMP  RATE  (V/s) 


(b) 


RAMP  RATE  (V/s) 


Figure  7.  Overvoltage  as  a  function  of  ramp  rate:  (a)  Calcu¬ 
lated  data  with  three  assumed  dc  breakdown  voltages; 
(b)  From  measured  characteristics  of  Harari  [11]. 


Figure  8.  Comparison  of  calculated  and  measured  over¬ 
voltage  as  a  function  of  ramp  rate.  Measured  curve  Is  for  a 
161-A  oxide,  whereas  calculated  curve  Is  for  500  A. 
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The  breakdown  current,  defined  as  the  current 
at  the  time  of  the  voltage  maximum,  was  also  found 
to  increase  with  ramp  rate.  This  is  also  observed 
for  Harari's  measurements,  shown  in  figure  2.  De¬ 
fining  breakdown  power  as  that  pertaining  to  the 
maximum  voltage  for  the  calculations  and  for  the 
breakpoint  in  Harari’s  measurements,  breakdown 
power  is  plotted  as  a  function  of  time  to  breakdown 
in  figure  9.  It  is  seen  that  breakdown  power  varies 
inversely  with  breakdown  time  for  the  ca'culations 
and  to  the  -0.9  power  for  the  measurements. 
There  is  an  order-of-magnitude  error  of  4  in  the 
powers  extrapolated  to  a  ramp  rate  of  1  V/s,  as 
shown  in  figure  10,  but  such  an  extrapolation  of  10 
orders  of  magnitude  is  probably  meaningless.  In  sil¬ 
icon,  the  slope  of  the  power  versus  breakdown  time 
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varies  from  -1  to  -1/2  and  less  in  only  a  few  orders 
of  time  magnitude. 

5.  Sinusoidal  Applied  Voltages 

Calculations  have  also  been  made  for  break¬ 
down  with  one-half  period  of  a  sinusoidal  wave¬ 
form.  Contrary  to  the  previous  practice  [12]  of  ad¬ 
ding  a  dc  bias  to  the  sinusoidal  voltage,  for  these 
calculations  the  computer  initial  conditions  were 
chosen  with  an  initial  time  greater  than  zero,  thus, 
an  applied  voltage  greater  than  zero.  This  assures 
an  accurate  sinusoidal  waveform.  Calculations 
were  made  for  the  100-nm  oxide  for  frequencies  of 
1.0  to  40  GHz.  Figure  1 1  shows  current  buildup  as 
a  function  of  time  for  applied  maximum  voltages 
just  below  and  above  the  breakdown  values.  It  is 
seen  that  the  breakdown  voltage  and  current  in¬ 
crease  with  frequency.  Figure  12  shows  the  lowest 
applied  power  to  cause  breakdown  during  the  first 
half-cycle  as  a  function  of  frequency.  Also  shown  is 
the  highest  applied  power  for  which  breakdown  did 
not  occur.  If  the  dissipated  power  in  the  diode  is 
plotted  against  time  to  breakdown,  as  was  done  in 
figure  9,  a  slope  of  -1.0  is  again  obtained,  and  the 
intercept  is  about  a  factor  of  three  lower  in  power 
for  the  sinusoidal  case. 

Earlier  sinusoidal  waveform  calculations  were 
made  with  a  large  dc  approach  voltage,  plus  a 
sinusoidal  increment.  The  effective  frequency  must 
be  estimated  by  fitting  a  true  sinusoidal  wave  to  the 


Log  Time  (s) 


Figure  9.  Breakdown  power  as  a  function  of  time  to  break¬ 
down:  (a)  calculated  for  a  100-nm  film  thickness  and 
(b)  from  measured  characteristics  of  Hararl  [11]  for  a 
16.1-nm  film. 


Figure  10.  Comparison  of  calculated  and  measured  break¬ 
down  power  as  a  function  of  breakdown  time.  Measured 
curve  Is  for  a  1 61  -A  oxide,  whereas  calculated  curve  is  for 
1000  A. 


Figure  11.  Current  density  as  a  function  of  time  for  various 
frequencies  as  labeled.  Applied  voltages  resulting  In  break¬ 
down  and  nonbreakdown  are  given  for  each  frequency. 


incremental  portion.  For  example,  for  an  applied 
voltage,  Va,  given  by 

Va«  180 +  80  sin  1  x  101  ’t  (2) 

a  good  fit  is  obtained  for  a  frequency  of  8.33  GHz, 
as  compared  to  15.9  GHz  for  the  sinusoidal  incre¬ 
ment.  For  the  calculation  made  with  Va  given  by 
equation  (2),  the  oxide  failed  to  break  down,  but 
with  an  85-V  sinusoidal  increment,  breakdown 
occurred.  A  second  half  cycle  was  calculated  with 
the  same  polarity  and  the  Va  of  equation  (2),  but 
with  the  distribution  of  holes  remaining  in  the  oxide 
at  the  end  of  the  first  half  cycle  (actually,  the 
distribution  was  obtained  at  a  voltage  somewhat 
above  the  approach  voltage  of  180  V).  Breakdown 
occurred  in  the  second  half  cycle  at  a  voltage 
about  4  V  below  that  for  breakdown  in  the  first  half 
cycle.  The  second  half-cycle  calculation  was 
repeated  with  the  opposite  polarity.  Breakdown  oc¬ 
curred  at  80  V  below  that  of  the  first  half  cycle.  Fur¬ 
ther  investigation  showed  that  ionization  by  the 
stored  holes  was  responsilbe  for  the  lowered 
breakdown  voltage.  Two  first  half-cycle  calcula¬ 
tions  were  repeated  with  the  hole  ionization  coeffi¬ 
cient  equal  to  zero,  and  the  breakdown  voltage  was 
virtually  unchanged. 
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Figure  12.  Calculated  applied  power  for  breakdown  as  a 
function  of  frequency.  Breakdown  Is  denoted  by  symbol  x 
and  nonbreakdown  by  symbol  o. 

It  is  highly  probable  that  ionization  by  holes  does 
not  occur,  but  since  the  transport  of  holes  is  not  un¬ 
derstood,  it  is  possible  that  ionization  could  occur. 
If  an  experimental  test  showed  a  lower  breakdown 
with  a  reversed  polarity,  ionization  by  holes  could 
be  indicated. 

6.  Discussion  and  Conclusions 

Calculations  have  been  made  of  breakdown  in 
Si02  films  with  the  assumption  of  avalanche  and 
recombination  processes.  Applied  voltages  were 
(1)  constant,  (2)  a  constant  ramp  rate,  or  (3)  sinu¬ 
soidal  in  form.  The  time  scale  of  the  calculations  is 
limited  to  about  100  ps.  There  is  a  10-orders-of- 
magnitude  difference  in  the  time  scales  used  in  the 
calculations  and  in  the  most  informative  exper¬ 
imental  measurements.  Despite  this  large  gap 
there  is  a  real  sense  of  agreement,  both  in  terms  of 
overvoltage  as  a  function  of  ramp  rate  and  in  power 
to  breakdown  as  a  function  of  time  to  breakdown, 
between  the  calculations  and  the  measurements. 
Others  have  reported  general  agreement  between 
bulk  laser  breakdown  fields  and  those  measure¬ 
ments  obtained  near  dc  in  oxide  films.  It  would  be 
very  helpful  if  measurements  could  be  made  on  an 
intermediate  time  scale. 
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